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THE  EULERIAN  TIME  CORRELATION  FUNCTION  IN 
HOMOGENEOUS  ISOTROPIC  TURBULENCE 

R.  RUBINSTEIN*  AND  GUOWEI  HEt 

Abstract.  Two  general  models  are  proposed  for  the  Eulerian  time  correlation  function  in  homogeneous 
isotropic  turbulence.  The  first  is  based  on  continued  fraction  approximations  to  its  Laplace  transform, 
and  the  second  is  based  on  random  sweeping  by  a  possibly  non-Gaussian  velocity  field.  Both  models  can 
give  reasonable  quantitative  agreement  with  DNS  data  for  moderate  time  separations  over  which  the  time 
correlation  functions  at  different  wavenumbers  exhibit  a  common  self-similar  form. 

Key  words,  time  correlation,  turbulence,  sub-Gaussian 

Subject  classification.  Fluid  Mechanics 

1.  Introduction.  The  analysis  of  time  correlations  in  turbulence  begins  with  Kraichnan’s  explanation 
[1]  of  the  inconsistency  of  Eulerian  turbulence  closures  with  Kolmogorov  scaling.  This  work  revealed  that 
the  dynamic  decorrelation  mechanisms  are  distinct  for  Eulerian  and  Lagrangian  time  correlations  and  that 
energy  transfer  in  turbulence  must  be  analyzed  in  terms  of  Lagrangian  quantities. 

Since  Lagrangian  time  correlations  arise  both  in  the  problem  of  energy  transfer  in  turbulence  and  in  the 
equally  fundamental  problem  of  passive  scalar  diffusion,  they  would  seem  to  be  more  important.  Nevertheless, 
Eulerian  properties  are  relevant  in  a  broad  class  of  problems  in  which  turbulence  acts  as  a  time-dependent 
random  medium  and  properties  at  fixed  locations  in  space  are  required.  Wave  scattering  by  turbulence  is 
one  example;  provided  the  turbulence  is  not  simply  frozen  during  the  passage  of  the  wave,  analysis  of  this 
problem  will  depend  on  Eulerian  time  correlations. 

Sound  radiation  by  turbulence  is  another  problem  in  which  Eulerian  correlations  might  be  relevant.  As 
applied  to  homogeneous,  isotropic  turbulence,  Lighthill’s  theory  [2]  appears  to  treat  fixed  regions  of  space, 
not  moving  volumes  of  fluid,  as  sound  sources.  Since  the  space-time  properties  of  fixed  spatial  volumes 
determine  the  sound  source,  Eulerian  time  correlations  are  relevant.  Since  the  observer  defines  a  coordinate 
system  at  rest,  Eulerian  time  correlations  are  certainly  consistent  with  Lighthill’s  theory.  More  generally, 
Kaneda  has  suggested  that  whereas  energy  transfer  in  turbulence  depends  on  Lagrangian  time  correlations, 
momentum  transfer  depends  on  Eulerian  time  correlations,  which  also  suggests  the  relevance  of  Eulerian 
time  correlations  to  sound  radiation. 

The  Eulerian  time  correlation  function  is  defined  in  homogeneous,  isotropic  turbulence  by 

(1.1)  (u(x,  t)  ■u{Tc  +  r,t  +  r))  =  Q{r)  tpir,  r) 
where 

(1.2)  (3(r)  =  (u(x,t)  •  u(x -|-r,t)) 
is  the  single-time  correlation  function. 
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The  basic  fact  about  Eulerian  time  correlations  is  the  similarity  form  for  inertial  range  separations  [1,  3] 


(1.3)  ip(r,T)  =  Ip 
where 

(1.4)  =  (u(x,t)  •  u(x,t)) 
is  the  rms  velocity  fluctuation. 


Fig.  1.1.  Time  correlation  function  at  various  wavenumbers  plotted  as  functions  of  the  similarity  variable  ^  =  Vkr.  The 
thick  curve  with  the  lowest  minimum  corresponds  to  a  low  wavenumber;  increasing  wavenumber  corresponds  to  increasing  ip{^) 
at  ^  =  4. 

The  validity  of  the  similarity  form  Eq.  (1.3)  has  been  demonstrated  repeatedly  in  the  literature  [3].  We 
confirm  this  similarity  form  with  new  observation  summarized  in  Eig.  1.1,  which  shows  the  time  correlation 
function  in  Eourier  space,  'tp{k,T)  for  various  wavenumbers  k,  plotted  as  functions  of  the  similarity  variable 
Vkr.  The  data  is  taken  from  He  et  al.  [4].  The  time  correlations  for  all  wavenumbers  collapse  very  well  for 
time  separations  Vkr  <  3.  At  longer  times,  the  time  correlation  functions  cross  zero  and  oscillate.  These 
features  appear  to  be  quite  robust,  although  our  data  is  not  sufficiently  well  resolved  to  draw  definitive 
conclusions  at  very  long  time  separations.  This  behavior  is  somewhat  in  contrast  to  the  data  of  [5],  however, 
we  provisionally  assume  it  is  correct.  Although  the  correlation  functions  at  different  wavenumbers  no  longer 
collapse  so  well  when  Vkr  >  3,  the  time  correlation  functions  for  wavenumbers  in  the  range  10  <  A:  <  35, 
appear  to  exhibit  a  common  similarity  form  even  in  the  neighbourhood  of  the  zero-crossing.  We  tentatively 
identify  these  wavenumbers  with  an  inertial  range,  although  the  relatively  low-resolution  (128^)  DNS  does 
not  contain  an  extended  region  of  unambiguous  Kolmogorov  scaling. 

In  this  paper,  we  explore  some  possible  analytical  forms  for  the  function  tp.  Apart  from  its  purely 
theoretical  interest,  there  are  applications  in  which  the  precise  analytical  form  of  the  time  correlation  is 
important.  In  astrophysical  applications  [6],  the  total  acoustic  radiation  from  stars  depends  surprisingly 
sensitively  on  the  assumptions  made  about  the  time  correlation  function.  A  similar  sensitivity  to  the  time 
correlation  function  was  found  by  Bertoglio  et  al.  [7]. 
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2.  Continued  fraction  methods.  Kaneda  et  al.  [5]  curvefit  turbulent  time  correlation  functions  by 
rational  functions,  using  the  Taylor  series  approximations  and  the  Fade  table.  Thus,  if  for  fixed  k,  we  write 

(2.1)  i>{k,T)  =  uJo{k)  +  ^a;2(A:)r^  +  ^^uJi{k)T^  H - 

then  the  Taylor  coefficients  are  evaluated  as 


(2.2) 


LOo{k)  =  tp{k,0)  =  1 

uJ2{k)  =  ^(^>0)  =  - 

^4(A:)  =  ^{k,0)  =  + 


(u(k,t)  ■  u(-k,t)) 
(u(k,t)u(k,t)) 
(ii(k,t)ii(-k,t)) 
(u(k,t)u(k,t)) 


The  time  derivatives  on  the  right  side  can  be  replaced  by  single-time  moments  of  the  velocity  field  by  applying 
the  Navier-Stokes  equations.  The  time  derivative  of  order  n  then  is  given  in  terms  of  a  single-time  moment 
of  order  2n.  This  fact  connects  the  time  correlation  function  to  intermittency  [8]. 

This  method  of  calculation  leads  to  very  good  results.  Here,  a  different  method  is  followed,  which  leads 
to  a  rational  approximation  of  the  Laplace  transform  of  the  correlation  function  [9] .  In  molecular  hydrody¬ 
namics,  the  correlation  function  is  expressed  as  a  continued  fraction  through  the  Zwanzig-Mori  projection 
operator  formalism  [10].  Recently,  this  formalism  has  been  generalized  to  non-equilibrium  conditions  by  [8]. 

However,  we  can  also  give  an  elementary  account  as  follows.  Define  the  Laplace  transform  of  the  time 
correlation  function  as  usual  by 

rOO 

(2.3)  ^(s)  =  /  dr  e“^’^^(r) 

Jo 

and  assume  the  continued  fraction  representation 


(2.4) 


V'(s) 


do  02  O4  Cl6 
S-b  S-b  S-b  S  + 


and  the  Taylor  series  expansion  Eq.  (2.1).  Substituting  Eq.  (2.1)  in  Eq.  (2.3), 


(2.5) 


,  .  ,  1  1-O2  1-O4 

ip{s)  =  -  -b  ^  +  -r  + 


Equating  the  expressions  for  ^(s)  Eqs.  (2.4)  and  (2.5),  the  Euclidean  algorithm  gives 


02  —  —LO2 

(2.6)  04  = - b  t02 

CO2 

CJ4  —  CJe  1102  1^4 

0-6  = - - 1 - 

lu2  ~  W4/a;2  C02 


An  obvious  approximation  method  is  finite  truncation  of  this  continued  fraction  [10].  We  obtain  in  this  way 
first  the  the  one-parameter  or  Markovian  model 


(2.7) 


V'(t)  =  e 


and  at  the  next  order,  the  two-parameter  model 


(2.8) 


ip{r)  =  ^cosw^  -b  ^  sinw^j 


e-«/2 
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In  Eqs.  (2.7)  and  (2.8),  ^  =  Vkr,  the  similarity  variable  corresponding  to  Eq.  (1.3) 

It  should  be  stressed  that  an  assumption  of  complete  similarity  is  made  in  writing  Eqs.  (2.7)  and  (2.8), 
namely  that 

(2.9)  L02n  ~ 

This  is  an  assumption  of  ‘normal  scaling’  for  the  Taylor  coefficients;  in  view  of  Eqs.  (2.6),  it  also  implies 
normal  scaling  of  the  coefficients  in  the  continued  fraction  expansion  Eq.  (2.4).  This  assumption  may  be 
satisfactory  for  Eourier  coefficients  [11]. 

The  Markovian  model  is  well-known  to  be  incorrect,  because  it  does  not  have  zero  slope  at  r  =  0.  But 
Eig.  1.1  also  suggests  that  the  two-parameter  model  is  qualitatively  inadequate,  because  it  predicts  equally 
spaced  zeroes  of  the  correlation  function.  We  could  proceed  to  the  general  three-parameter  model,  but  prefer 
to  consider  a  special  form  of  this  model  [10]  defined  in  terms  of  the  memory  function 


(2.10) 

K{t)  =u;‘^{l  +  LUT)e-^^ 

by  the  equation 

(2.11) 

+  f  dr' K(t  —  =  0 

JO 

Eq.  (2.11)  is  equivalent  to  a 

forms  as 

system  of  ordinary  differential  equations,  expressed  in  terms  of  Laplace  trans- 

(2.12) 

S1p{s)  +  =  1 

s^l(s)  -1-  LO(l>i{s)  -  'tp{s)  -  LO<j)2{s)  =  0 

S^2(s)  +  i^(t>2{s)  -  tp{s)  =  0 

which  proves  more  convenient  for  numerical  evaluation. 

Since  this  is  a  two-parameter  model,  it  is  defined  by  two  properties  of  the  correlation  function,  which 

we  take  to  be  the  microscale 

(2.13) 

and  the  integral  scale 

^2  1 
"  m 

(2.14) 

POO 

tm=  dt  ip(t) 

Jo 

The  connection  between  these  quantities  and  the  model  parameters  is 

(2.15) 

(2.16) 

~2  ^ 

LU  =  - 

u;  tm 

The  generality  of  this  class  of  functions  is  evident  from  the  equivalent  system  of  equations  Eq.  (2.12), 
which  shows  that  pure  exponential  decay  occurs  in  the  limit  w  ss  0  and  that  oscillations  occur  in  the  limit 
a;  PS  0.  Intermediate  behavior  can  be  anticipated  between  these  limits  and  is  confirmed  in  Eig.  2.1.  Trial  and 
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Fig.  2.1.  The  two-parameter  eorrelation  funetion  of  Eqs.  (2.10)-(2.11)  plotted  for  various  values  of  the  ratio  r  =  uj/oj: 
r  =  0  (solid),  r  =  1  (dashed),  r  =  10  (dot-dash).  As  r  inereases  from  zero,  the  funetion  goes  from  pure  deeay  to  pure  oseillation. 


Fig.  2.2.  The  eorrelation  funetion  Eqs.  (2.10)-(2.11)  fit  to  the  data  of  Eig.  1.1.  Symbols  are  DNS,  the  solid  line  is  the 
eurvefit,  and  the  dashed  line  is  the  best  squared  exponential  eurvefit.  The  fitting  was  done  by  trial  and  error.  The  region  of 
negative  values  of  the  eorrelation  funetions  is  reprodueed  well  by  the  empirieal  funetion. 


error  adjustment  of  the  parameters  led  to  an  approximation  to  the  measured  Eulerian  correlation  functions 
which  is  shown  in  Fig.  2.2.  We  emphasize  that  the  relations  Eqs.  (2.15)-(2.16)  were  not  used  to  obtain  the 
eurvefit. 

This  result  shows  that  a  good  fit  to  the  time  correlation  function  can  be  obtained  in  principle  knowing 
only  the  microscale  and  integral  scale.  Nevertheless,  the  fact  remains  that  neither  the  equations  relating  the 
Taylor  coefficients  and  the  model  parameters  Eq.  (2.6)  nor  the  system  of  ordinary  differential  equations  for 
the  correlation  function  Eq.  (2.12)  has  much  apparent  physical  significance. 

3.  The  random  sweeping  model.  An  important  theoretical  model  for  the  Eulerian  time  correlation 
function  is  provided  by  the  direct  interaction  approximation  [12].  However,  the  functional  form  predicted, 
V'(t)  ~  JiiOK  i®  unrealistically  oscillatory. 

Another  approximation  for  Eulerian  time  correlations  is  based  on  Kraichnan’s  model  problem  for  the 
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random  sweeping  effect  [1].  In  this  problem,  a  fixed  frozen  random  velocity  field  u(x)  is  advected  by  a 
random  constant  velocity  field  v.  The  result  is  the  time  dependent  field  defined  by 

(3.1)  u(k,  t)  =  — fk  •  vu(k,  t) 
where 

(3.2)  u(k,0)=u(k) 
so  that 

(3.3)  u(k,t)  =  exp(iv  •  kt)u(k) 

If  V  is  Gaussian,  then 

(3.4)  ^(r)  =  (exp(iv  •  kr))  =  exp(— G^r^A:^/2) 
where 


(3.5) 


=  (v  •  v) 


The  time  correlations  of  the  actual  turbulent  velocity  field  are  approximated  by  the  (ensemble-averaged) 
time  correlations  of  the  velocity  field  defined  in  Eq.  (3.3)  The  assumption  underlying  this  model  is  that  the 
large  scales  responsible  for  sweeping  in  the  actual  velocity  field  can  be  considered  independent  of  the  small 
scales  at  which  the  time  correlation  is  evaluated.  The  first  equality  in  Eq.  (3.4)  shows  that  the  pdf  is  related 
to  the  frequency  spectrum. 

It  is  evident  that  this  model  will  not  provide  a  very  good  fit  to  the  data  because  the  measured  correlation 
functions  all  take  negative  values.  We  consider  how  this  model  might  be  generalized  to  obtain  better 
agreement.  If  the  pdf  of  v  is  non-Gaussian,  then  certainly  the  second  equality  in  Eq.  (3.4)  does  not  apply. 
There  is  considerable  recent  work  [13,  14]  suggesting  that  the  single-point  pdf  of  velocity  fluctuations  is  not 
exactly  Gaussian,  but  is  instead  sub- Gaussian:  that  is,  the  tails  are  lighter  than  Gaussian,  in  contrast  to 
familiar  ‘intermittency’  behavior. 

A  simple  sub-Gaussian  is  defined  by  truncating  a  Gaussian, 


(3.6) 


p(v) 


6exp(— u^/2E^)  if  u  <  a 
0  if  u  >  a 


where  the  standard  deviation  oc  V  defines  the  relevant  sweeping  velocity.  Limiting  cases  of  the  corresponding 
correlation  function  are  the  squared  exponential  Eq.  (3.4)  when  a  =  oo  and 

(3.7)  ^ 

when  a  =  0.  In  general,  this  correlation  function  satisfies  the  inhomogeneous  differential  equation 


(3.8) 


dilf 

dr 


+  =  —TTO^e 


d  sin(ar) 
dr  ar 


Eig.  3.1  plots  this  correlation  function  for  various  values  of  the  parameter  a. 

Although  these  correlation  functions  look  promising,  there  is  little  doubt  that  the  underlying  model  for 
the  pdf  is  somewhat  artificial.  Ealkovich  and  Lebedev  [13]  have  attempted  to  compute  the  tail  of  the  single¬ 
point  velocity  pdf.  The  tail  is  found  to  depend  on  the  forcing  statistics,  and  is  therefore  not  universal.  Simply 
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to  illustrate  some  possibilities,  we  follow  this  article  and  consider  a  non-Gaussian  pdf  with  an  exp(— tail. 
Suppose,  to  have  consistency  with  near  Gaussian  behavior  for  small  deviations  from  zero,  that 


(3.9) 


p(v)  =  aexp[— +  av‘^)l2V‘^] 


Then  satisfies 


(3.10) 


(fip  2  dtp 
T  dr 


1  dip 

2  dr 


1 


Tip  =  0 


Note  that  when  a  =  0,  the  time  correlation  function  reduces  correctly  to  the  squared  exponential  form 
corresponding  to  a  Gaussian-distributed  sweeping  velocity.  Two  typical  forms  of  this  function  are  shown  in 
Fig.  (3.2). 


Fig.  3.1.  The  one-parameter  eorrelation  funetion  defined  by  Eq.  (3.8)  plotted  for  various  values 
a  =  100  (solid),  a  =  10  (dashed),  a  =  1  (dot-dash).  The  limiting  eases  are  =  exp(— when  a  =  oo 


of  the  parameter  a: 
and  =  sin(G/^ 


when  a  =  0. 


We  note  the  existence  of  asymptotic  solutions  to  Eq.  (3.10) 

(3.11)  tpir)  ~  cos('s/3^"‘/^/2  +  <p)  exp(— ^"‘/^/2) 

This  solution  indicates  that  the  correlation  function  falls  off  considerably  more  slowly  than  exp(— ^^)  and 
that  the  correlation  function  crosses  zero  with  increasing  frequency  as  r  oo.  In  this  respect,  the  simpler 
model  Eq.  (3.6)  seems  more  realistic,  but  the  fact  is  that  the  available  data  do  not  rule  Eq.  (3.11)  out. 

In  Eig.  (3.3),  this  function  is  compared  to  the  empirical  distribution  corresponding  to  a  variety  of 
wavenumbers  k.  Again,  the  curvefitting  was  done  by  trial  and  error,  with  no  attempt  to  determine  the  best 
fit  from  numerical  properties  of  the  measured  correlation  function.  A  best  fit  squared  exponential  is  included 
for  comparison.  This  figure  shows  that  the  low  k  correlation  function  is  fit  remarkably  well  by  this  type  of 
function. 

Despite  the  good  fit  for  the  low  wavenumber  case,  these  sub-Gaussian  models  are  generally  too  ‘elastic’ 
in  comparison  to  the  data:  the  negative  amplitude  after  the  first  zero  crossing  tends  to  be  larger  than  the 
data,  and  it  tends  to  recover  too  quickly.  This  raises  the  question  of  the  accuracy  of  the  data  in  this  range. 
More  accurate  measurements  should  be  made  in  order  to  make  more  definitive  comparisons. 
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0  10  20 
time  separation 


Fig.  3.2.  The  one-parameter  eorrelation  funetion  defined  by  Eq.  (3.10)  plotted  for  two  values  of  the  parameter  a:  a  =  100 
(solid),  a  =  10  (dotted). 


time  separation  time  separation 

Fig.  3.3.  The  eorrelation  funetion  Eq.  (3.10)  fit  to  the  data  of  Eig.  1.1.  Symbols  are  DNS,  the  solid  line  is  the  eurvefit, 
and  the  dashed  line  is  the  best  Gaussian  (squared  exponential)  eurvefit.  The  graph  on  the  left  eorresponds  to  relatively  high 
wavenumbers,  the  graph  on  the  right  to  a  low  wavenumber.  The  Gaussian  provides  a  good  fit  to  the  high  wavenumber  data, 
but  at  low  wavenumber,  the  sub- Gaussian  provides  mueh  better  agreement. 

4.  Conclusions.  The  collapse  of  the  Eulerian  time  correlation  function  to  a  self-similar  form  based  on 
the  sweeping  velocity  is  confirmed  by  DNS  data.  For  moderate  time-separations,  the  collapse  is  excellent 
for  all  wavenumbers,  and  continues  for  larger  time  separations  for  inertial  range  modes.  Models  for  the  time 
correlation  function  based  on  continued  fraction  approximation  and  on  a  random  sweeping  model  can  both 
give  satisfactory  agreement  with  the  data.  More  accurate  measurements  of  the  time  correlations  at  large 
time  separation  will  be  required  to  refine  these  models. 
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